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THE  STATEMENT  Of  THK  PROBLEM  OT  THiS  STATIONARY  .STRUCTURE  OF  THE 
BOITODARY  LAYER  OF  THE  A-IMOSUHERE 

by 

1).  li,  Lalkhtmnn 

Th*  overall  problam  of  the  distribution  of  wind,  teiaonrature, 
humidity  and  of  the  turbulence  coefficient  in  the  boundary  layer 
of  the  atmosphere  under  stationary  conditions  is  mathemattcfilly 
formulated  in  thia  article. 

The  rasulta  of  the  calculations  are  compared  with  experimental 

data. 

When  an  air  masn  moves  above  a  homogeneous  underlying  surface  as 
a  result  of  turbulent  mixing,  the  determined  profiles  of  temiiernturo , 
humidity,  wind  velocity,  and  of  the  turbulence  coefficient  are  ostab- 
llahad.  In  a  stationary  condition  all  the  charaoterlsltcs  mentioned 
are  Interrelated,  This  connection  la  caused  first  of  all  by  processae 
on  the  active  surface,  as  a  result  of  which  a  fixed  quantity  of  radiant 
energy  is  distributed  along  three  channels:  to  the  warming  of  the  noil 
and  air,  and  to  evaporation  and  secondly  by  the  machaniam  of  turbulence 
itself  which  ia  characterized  by  the  fact  that  the  intenaity  of  turbulent 
mixing  and  the  vertical  temperature  gradients  and  the  wind  velocities 
condition  each  other,  In  so  far  as  the  turbulent  energy  influx  due  to 
the  energy  of  average  motion  and  its  expenditure  or  work  agalnat  Arch- 
Imedlan  forces  sre  determined  by  the  vertical  gridlonta  of  temperature 
and  wind  velocity.  At  the  some  time,  the  Intenni flcation  or  weakening 
of  turbulent  mixing  exerts  a  rapid  and  a  completely  determined  effect  on 
the  vertical  profile  of  any  substance.  In  this  connection,  the  correct 
slateraent  of  the  problan  of  the  etatlonary  metooroiogl cni  regime  of  the 
boundary  layer  of  the  atmosphere  muet  be  the  common  solution  of  a  certain 
system  of  aquations  for  determining  the  t  empoi'stu  re ,  humidity,  wind 
velocity,  and  turbulence  coeTficient, 
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Only  four  equations  ore  known  to  find  the  five  unknown  functions j 
T,  u,  T,  k,  qt  where  T  ia  the  air  temperature,  q  is  the  specific  hu»i- 
dlty  of  the  air,  u,  w,  ajr-e  the  components  of  the  wind  velocity,  k 
ia  the  turbulence  coefficient;  they  sres  two  equations  of  motion,  the 
equation  of  heat  conductivity,  and  the  aquation  for  the  diffualon  of 
water  vapor.  In  thia  connection,  it  was  necessary  to  carry  out  a 
significant  number  of  investigationa  with  the  assumption  that  the 
coefficient  of  turbulence  is  a  known  function  of  height. 

Although  important  results  of  a  pui'ely  practical  aspect  were 
obtained  in  a  series  of  works,  they  yielded  comparatively  little  for 
theory.  An  explanation  of  the  variations  in  the  vertical  profiles  of 
meteorological  elements  by  any  changes  in  the  intensity  of  turbulent 
exchange,  which  result  from  such  investigations,  ia  hardly  satisfac¬ 
tory.  The  phenomena  mentioned  are  connected  not  as  causa  and  effect, 
but  rather  aro  the  manifestation  of  a  single  process,  and,  only  because 
of  this,  are  related  to  each  other.  On  the  same  basis,  It  is  possible 
to  explain  the  variations  of  the  temperature  profile  by  the  change  of 
turbulence,  and  to  explain  changes  in  the  intensity  of  turbulent 
exchange  by  changes  of  the  temperature  gradients.  A  considerably  more 
correct  statement  is  the  application  to  the  problem  of  methods  of  the 
theory  of  similarity  and  of  an  analysis  of  the  measuroments ,  which 
has  been  made  recently  in  a  aeries  of  investigations.  Under  the  present 
condition  of  the  problem,  this  method  ia  completely  logical  and  has 
already  made  It  possible  to  obtain  a  series  of  rcsultn  important 
for  understanding  the  mechanism  of  the  phenomenon.  Unfortunately, 
with  a  more  or  less  complete  statement  of  the  problem,  the  results  are 
indefinite 

The  application  of  the  energy  bolance  of  turbulence  na  a  non- 
satiafylng  equation  aeema  fully  expedient  to  us.  One  can  easily 
underotand  that  those  name  unknown  functions  enter  thia  equation 
and  it  la  not  the  result  of  the  four  mentioned  equations.  However 
the  application  of  the  indicated  equation  has  many  dl  f  f  i  c\il  1 1  on  .  First 
of  nil,  there  are  dlfficultlso  of  n  principle  order  that  nrlno  from  the 
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fact  that  expressions  for  the  diffusion  of  turbulence  energy  and 
for  its  dissipation  into  heat  are  presently  unknown.  Important  diff¬ 
iculties  of  a  mathematical  nature  also  arise  with  a  correct  state¬ 
ment  because  t)ie  system  of  equations  i a  the  indicated  statement 
becomes  essentially  uon-linear. 

Let  ns  now  indicate  a  certain  method  which  makes  it  possible  to 
obtain  a  number  of  important  results  in  the  solution  of  the  general 
problem.  Let  us  uae  Lho  fact  that  it  is  poseiblo  to  approximate  the 
turbulence  coafficioat  aa  a  function  of  height  as  follows: 


k{2) 


'  fcrr.s^/i, 
)'  ’  ^ 


Vftien  there  ia  equilibrium  of  stratification  c  ■  0,  during  the 
steady  state  c  >  0,  and  during  the  unsteady  state  c  <  0, 

First  let  us  examine  the  statement  of  the  problem  for  the 
equilibrium  state.  If  one  is  to  consider  that  for  sufficiently 
BBali  values  of 


^■(C,)  Ay 

In 

where  is  the  modulus  of  the  wind  velocity  at  height  is 

the  roughness  of  the  underlying  surface,  then  instead  of  the  unknown 
function  k(E),  only  the  height  of  the  cleavage  h  is  unknown.  In 
this  connection,  we  can  present  the  equation  of  the  bnlanco  of  ths 
energy  of  turbulence  in  Integral  form,  after  having  integrated  it 
with  respect  to  the  entire  boundary  layor.  In  such  a  cane,  we  can 
disregnrd,  bnenune  of  its  relative  amailnese,  thn  comj)on»nt  which 
takoa  into  account  the  diffusion  of  Lho  energy  of  turbulence. 
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ka  to  the  dlaaipation  of  the  turbulence  energy  Into  heat, 
we  can  obtain  a  oonveninnt  formula  for  it  by  taking  advantage  of 
certain  phyaical  considemtiona  of  the  dioaipation  raeohaniaa,. 

The  dissipation  of  turbulence  energy  into  heat  ia  equal  to  the 
work  of  the  forcea  of  friction,  which  occurs  because  the  elements 
of  turbulence  move  with  a  velocity  that  ia  different  than  the 
velocity  of  the  surrounding  medium. 

We  shall  introduce  ths  following  daolgnationa i  T  is  the  force 
of  friction,  p  ia  the  air  density,  w*  is  the  velocity  pulsation, 
characteristic  for  a  certain  vortex,  v  Is  the  kinenatlo  viscosity, 

£  ia  the  characteristic  dimension  of  the  vortex. 

Obviously  /,  w'.  V)  ---  ■ 


The  work  of  the  force  of  friction  per  unit  time  is  fw',  and 
the  work  attributed  to  a  unit  of  volume  and  equal  to  the  dissipation 
of  the  turbulence  energy  into  heat,  can  be  expressed  par  unit  of 
tine  per  unit  of  volume  as 


(1) 


All  values  entering  into  (1)  are  averaged.  If  we  take  advantage 
of  the  fact  that 


and 


r  dt 


k  -  /  w' 


then 


(2) 
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In  such  a  caatf  the  balance  of  the  turbulence  energy  can  ba 
written  in  the  following  fora: 


dt 


,  dv 
\d: 


M. 


(3) 


M  is  the  diffueion  of  turbulence  energy. 

Here  u,  t  are  the  two  components  of  the  vector  of  the  averaged 
velocity. 

Having  integrated  aquation  (j)  with  respect  to  the  entire 
boundary  layer,  for  the  given  condition  we  will  obtains 
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(fz  ■==  0  —  M. 


(4) 


Thus  the  problem  of  finding  the  average  profiles  of  the  meteoro¬ 
logical  elements  and  of  the  coefficient  of  turbulence  in  the  boundary 
layer  of  the  atmosphere  reduces  to  the  integration  of  the  following 
system  of  equations: 
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Here  8  ia  the  aoil  tenparature,  ^  la  the  depth  measured  from 

the  surface  of  the  aollj (  1  >  6  )  Is  the  averaged  value  of  the 

dissipation  factor,  H  is  the  height  of  the  boundary  layer. 

If  w«  are  to  consider  that  the  radiant  flux  la  equal  to  zero 

or  is  given  in  the  form  of  some  interpolated  formula,  and 


I  (-*  Y  ‘  for  z  n 

liir 


(13) 


then  to  find  the  eight  unknowns,  u,  v,  T,  q,  k,  H,  b,  6  we  obtain 
eight  squationsj  the  differential  equations  being  linear. 

'Ctit  boundary  conditions  are  completely  clear;  on  the  surface 
layer.  It  is  necessary  to  fulfill  the  equation  of  the  heat  balance, 
and  to  give  the  temparMture  continuity  and  the  ela.^ticity  of  the 
water  vapor;  at  the  upper  edge  of  the  boundary  layer  It  is  natural 
to  give  the  tenperfiture ,  the  humif'lty  end  value  of  the  geostrophlc 
wind,  and  on  the  ground  it  la  neceseary  to  give  the  temperature  at 
the  depths  where  the  annual  variation  practically  disappears.  These 
boundary  conditions  provide  «  well-defined  aolution  for  the  problem. 

Wc  can  solve  the  formulated  problem  for  a  more  particular,  but 
auffioicntly  Important  case.  The  queatlon  is  the  distribution  of 
temperature,  wind  velocity,  and  of  the  turbulence  coefficient  uader 
different  conditions  of  vertical  stability.  Then  the  system  is  closed 
by  equations  (?),  (6),  (7),  (8),  (12)  and  (13).  In  auch  a  atateaeni, 
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th«  prob}.em  ia  •xauninad  in  datail  in  the  work  of  L.  B.  Orlanko 
and  Q.  Kh.  Taeitln  (  Tasitlc,  Q.  Kh«  and  L.  R,  Orlenko,  "Stat- 
alo&srnoe  raapradalenie  iratra,  tamperatwry  1  turbulentnogo  obaena 
r  pogranichnom  aloe  voadukha  pri  rasliohnykh  aoatoianiiakh  uatolch- 
Irosti"  (Statlonar-y  diatribution  of  the  wind,  temperature  and  tur¬ 
bulent  exchange  in  the  boundary  layer  of  the  air  under  various 
stability  conditions),  QlaToaia  Qaoflgicheakala  Obeervatoriia  Trudy. 

Ho.  9^1  8-28,  i960  C  AHS-T-B-387-t-]  )•  The  diaaipation  factor  la 
detsrainod  on  the  basis  of  experimental  data  on  the  dependence  of 
Cj^/D  on  c  , 

On  the  basis  of  the  calculated  profile  kCn)  and  o(k)  we  can 
detersine  the  wind  strength  by  the  foreula 

Obviously,  the  average  period  of  the  gusts  at  different  heights 
can  be  evaluated  by  the  formula  t  ~  1  /(do/  daO 

figure  1  given  the  results  of  the  calculations  of  the  average 
dlstrlbntion  of  wind  velocity  and  direction,  and  of  the  oharacterie- 
tics  of  strength  and  of  the  turbulence  coefficient  which  Q.  Kh,  Tseltin 
and  L,  R,  Orlsnko  carried  out  according  to  the  scheme  cited  above 
for  roughness  1.8  om.  The  mean  monthly  raluee  of  the  wind  velocity 
and  direction  and  of  the  temperature  obtained  from  observations  In 
Pavlovak  ars  also  givsn  there.  The  good  agreement  between  experi¬ 
mental  and  observed  values  points  to  the  correctness  of  the  aeaumed 
physical  conceptions. 
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^4  a  -1  L-J  Hit  .  •  .  (  I  I  141  1  I  I  i  1  I  I  Li  i  u  iiJ  ^ 

^  9  /;  •  /,r  /5«  20^  py  0 

7.  4# 


(«)  -  wind  rolocity  (B/aec),  (b)  -  wind  direction  (deviation  from  the 
direction  of  the  wind  at  level  200  a),  (c)  air  temperature  (deviatlona 
from  the  surface  valuoa),  (aoUd  lines  are  the  moaaured  valueo  - 
dashed  lines  are  the  calculatetl  values);  (d)  is  tho  coefficient  of 
turbulence  (m/aeo),  (e)  is  tho  average  wind  pulsation  («/aec)|  1- 
is  the  steady  state,  2-  Is  the  unstable  state,  3  -  i"  the  equilibrium 
state  > 


Figure  1 
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